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Abstract. We study the existence of soliton and black hole solutions of four- 
dimensional su(N) Einstein- Yang-Mills theory with a negative cosmological constant. 
We prove the existence of non-trivial solutions for any integer N, with N — 1 gauge 
field degrees of freedom. In particular, we prove the existence of solutions in which 
all the gauge field functions have no zeros. For fixed values of the parameters (at the 
origin or event horizon, as applicable) defining the soliton or black hole solutions, if 
the magnitude of the cosmological constant is sufficiently large, then the gauge field 
functions all have no zeros. These latter solutions are of special interest because at 
least some of them will be linearly stable. 
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1. Introduction 

Since the discovery of non-trivial, spherically symmetric, soliton [lj and black hole 
[2] solutions of the su(2) Einstein- Yang-Mills (EYM) equations in four-dimensional 
asymptotically flat space-time, there has been an explosion of interest in the properties 
of solutions of EYM theory and its variants (see, for example, [3] for a review). The 
original soliton solutions were a surprise discovery because there are no gravitational 
solitons [1] nor solitons in pure Yang-Mills theory in flat space-time [5]. The black hole 
solutions are counter-examples to the "no-hair" conjecture [6], in that they have no 
global charge, and are indistinguishable from a standard Schwarzschild black hole at 
infinity. However, while the "letter" of the no-hair conjecture is violated, its "spirit" 
remains intact because these black hole solutions, like the solitons, are linearly unstable 

Introducing a negative cosmological constant changes the picture completely (see 
[8] for a recent review). In asymptotically flat space, soliton and black hole solutions 
exist at discrete points in the parameter space, but in four-dimensional asymptotically 
anti-de Sitter (adS) space, su(2) EYM solutions are found in continuous open sets of 
the parameter space [HI [Hll EH] • Of particular interest are those solitons and black holes 
for which the single gauge field function u has no zeros. In asymptotically flat space, 
all solutions are such that u must have at least one zero, and the number of unstable 
perturbation modes of the solutions is 2k, where k is the number of zeros of uo [12J. 
In asymptotically adS space, at least some (but, interestingly, not all pj]) solutions 
where u(r) has no zeros are linearly stable under spherically symmetric [9], [TQl [TT] and 
non-spherically symmetric [13] perturbations, provided |A| is sufficiently large. Stable 
black holes in this model violate the "no-hair" conjecture in that, at infinity, they 
are indistinguishable from magnetically-charged Reissner-Nordstrom-adS black holes. 
However, only one extra parameter (which can be taken to be the value of the gauge 
field function at the event horizon) is required to completely determine the geometry 
exterior to the event horizon. While this is not a global charge measurable at infinity, 
one might still argue that the "spirit" of the "no-hair" conjecture is still valid as only a 
small number of parameters are required to fully describe stable black holes. 

Recently soliton and black hole solutions of four-dimensional su(N) EYM theory 
with a negative cosmological constant have been found [H [HI [15] . For purely magnetic 
solutions, the gauge field is described by A^ — 1 gauge field functions Uj. As in the 
su(2) case, there are solutions for which all the gauge field functions have no zeros. 
For |A| sufficiently large, at least some of these are stable, for any value of A" [TH H6] . 
Stable black hole solutions in su(N) EYM in adS therefore require an additional N — 1 
parameters to completely describe the configuration, even though the metric is identical 
to magnetically charged Reissner-Nordstrom-adS at infinity. Therefore an unbounded 
number of parameters is required to determine the structure of stable hairy black holes 
in adS. 

Our purpose in the present paper is to prove analytically the existence of solitons 
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and hairy black holes in four-dimensional su(N) EYM in adS, confirming the numerical 
work presented in [15J. Existence theorems for the original soliton and black hole 
solutions of su(2) EYM in asymptotically flat space were challenging to prove [TFJ [18] . 
but it has been shown analytically in this case that there are an infinite number of 
soliton and black hole solutions of the field equations, parameterized by the radius of 
the event horizon (if there is one) and k, the number of zeros of the gauge field function 
u. It has also been proven that k > 0. For su(iV) EYM in asymptotically flat space, 
rather less analytic work exists in the literature [T9], [20] . Arguments for the existence of 
solutions for arbitrary A" are presented in [20] , but the only complete analytic work is for 
su(3) [19] . For su(2) EYM in adS, proving the existence of solutions |9] is considerably 
easier than in the asymptotically flat case, mostly due to the boundary conditions at 
infinity It is therefore reasonable to imagine that the corresponding proof for su(N) 
EYM in adS will also be significantly easier than in the asymptotically flat case. 

The outline of this paper is as follows. We start, in Section [21 by outlining the field 
equations, gauge field ansatz and boundary conditions for soliton and hairy black hole 
solutions of su(N) EYM theory with a negative cosmological constant in four space-time 
dimensions [13] . Some trivial solutions of the field equations are discussed in Section 
I2.3[ the most important of which (for our later analysis) are embedded su(2) solutions. 
We also review the salient features of the numerical solutions presented in [15] . The field 
equations, presented in Section 12.11 are singular at the origin, black hole event horizon 
(if there is one) and at infinity. The first part of the proof of the existence of solutions 
is therefore to prove local existence of solutions satisfying the appropriate boundary 
conditions (which are given in Section 12.21) at these singular points, and this is the 
subject of Section [3j It turns out that this is the most lengthy and technically intricate 
part of the whole existence proof, particularly for local solutions near the origin. Our 
method in Section [3] follows a well-established approach [TTt |20~| I2TI [22] which involves 
casting the field equations in an appropriate form, so that a standard theorem (Theorem 
ED [H]) can be applied. It is deriving a suitable form for the equations near the origin 
which is the most complicated part of the analysis. The standard theorem then gives us 
local existence of solutions in a neighbourhood of the singular point (i.e. the origin, event 
horizon or infinity), and an important consequence of this theorem is that the solutions 
are analytic in the parameters which determine them. This property of analyticity is 
central to the arguments, presented in Section HJ which prove the existence of genuinely 
su(N) soliton and black hole solutions of the field equations. 

Our existence theorems are presented in Sections 14.31 and 14.41 In Section 14.31 we 
show existence of su(N) EYM solutions in a neighbourhood of any embedded su(2) 
solution, and, as a corollary of this, the existence of su(N) solutions for which all the 
gauge field functions have no zeros for any value of |A|. The behaviour of the solution 
space for large |A| is discussed in Section 14.41 and, in particular, we show that for 
any fixed values of the "initial" parameters (defined in Section 14.41) determining the 
solution of the field equations, for sufficiently large |A|, the solution generated by these 
parameters will be such that all the gauge field functions have no zeros. Finally, our 
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conclusions are presented in Section 



2. su(iV) Einstein- Yang-Mills theory 

In this section we will describe in detail our field equations and boundary conditions 
for both soliton and black hole solutions. We will also describe some simple embedded 
solutions of the theory, as well as reviewing the key features of the numerical solutions 
studied in [151 . 



2.1. Ansatz and field equations 

We consider four-dimensional su(N) EYM theory with a negative cosmo logical constant, 
described by the following action, in suitable units: 

Seym = \ j d A x^ [R - 2A - Tr F^F^] , (1) 

where R is the Ricci scalar of the geometry and A the cosmological constant. Throughout 
this paper, the metric has signature (— , +, +, +) and we use units in which 4nG = 1 = c. 
In this article we focus on a negative cosmological constant, A < 0. Varying the action 
([1]) gives the field equations 

T^u Rfiu ~^Rq\iv A^^y, 

=D^ = V^/+[A M ,F/]; (2) 
where the YM stress-energy tensor is 

TV = Tr F M F„ A - ^,Tr F Xa F x °. (3) 

In equations ([IH3]) we have employed the usual Einstein summation convention, where a 
summation over repeated indices is understood. However, from now on all summations 
will be given explicitly, and this will be particularly important in Section [3J 

In this paper we are interested in static, spherically symmetric soliton and black hole 
solutions of the field equations (j2J), and we write the metric in standard Schwarzschild- 
like co-ordinates as: 

ds 2 = -liS 2 dt 2 + ^ dr 2 + r 2 d6 2 + r 2 sin 2 9 d<f) 2 , (4) 

where the metric functions \x and S depend on the radial co-ordinate r only. In the 
presence of a negative cosmological constant A < 0, we write the metric function \x as 

, , 2m(r) Ar 2 , . 

tir) = 1 " - — • (5) 

We emphasize that we are considering only spherically symmetric black holes and not 
topological black holes which have been found in the su(2) case [23J. 

The most general, spherically symmetric, ansatz for the su(N) gauge potential is 

A = Adt + Bdr + ^(C-C H ) dQ- % - [(C + C H ) sin6 + DcosO] d<f>,(6) 
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where A, B, C and D are all (N x N) matrices and C H is the Hermitian conjugate of 
C. The matrices A and B are purely imaginary, diagonal, traceless and depend only on 
the radial co-ordinate r. The matrix C is upper-triangular, with non-zero entries only 
immediately above the diagonal: 

<7 iJ41 = Wi (r)e*»« (7) 

for j — 1, . . . , N — 1. In addition, D is a constant matrix: 

D = Diag (N - 1, iV - 3, . . . , -N + 3, -N + 1) . (8) 

Here we are interested only in purely magnetic solutions, so we set ^4 = 0. We may also 
take B = by a choice of gauge [21]. From now on we will assume that all the Uj(r) are 
non-zero (see, for example, [25] for the possibilities in asymptotically flat space if this 
assumption does not hold). In this case one of the Yang-Mills equations becomes [21] 

7 . = Vj = l,...,7V-l. (9) 

Our ansatz for the Yang-Mills potential therefore reduces to 

A= l -{C -C H ) old - % - [(C + C H ) sin 6 + Dcosd] d<f>, (10) 
where the only non-zero entries of the matrix C are 

(ii) 

The gauge field is therefore described by the N — 1 functions 0Uj(r). We comment that 
our ansatz ([TO]) is by no means the only possible choice in su(N) EYM. Techniques for 
finding all spherically symmetric su( N) gauge potentials can be found in J2B] , where all 
irreducible models are explicitly listed for N <6. 

With the ansatz fllOp . there are N — 1 non-trivial Yang-Mills equations for the iV — 1 
gauge field functions ujf 

( 2Ar 3 \ 
rW' + [2m- 2r 3 p e — w\ + WjUj = (12) 



3 

for j — 1, . . . , N — 1, where a prime ' denotes d/dr 



I N 

-Lj2\M-^ l -N-l + 2 j y], (13) 



4r 4 
j'=i 

Wj=l- to* + \ {u]_, + u] +1 ) , (14) 
and c^o = — 0. The Einstein equations take the form 

m=fi& + rpe, — = — , (15) 

o r 



where 



N-l 
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Altogether, then, we have N + 1 ordinary differential equations for the AT + 1 unknown 
functions m(r), S(r) and ojj{r). The field equations (I12II15P are invariant under the 
transformation 

ujj(r) -> -ujj{r) (17) 
for each j independently, and also under the substitution: 

j^N-j. (18) 

2.2. Boundary conditions 

We are interested in black hole and soliton solutions of the field equations (j!2)15p . 
However, the field equations are singular at the origin r = if we have a regular, soliton 
solution, at an event horizon r = if there is one, and at infinity r — > oo. Since the 
cosmological constant A < 0, there is no cosmological horizon. We need to derive the 
boundary conditions at each of these singular points. For A = 0, local existence of 
solutions of the field equations near these singular points has been rigorously proved 
[22]. The extension of these results to A < will be the focus of Section [31 



2.2.1. Origin The boundary conditions at the origin are the most complicated of the 
three singular points. We postpone the detailed derivation of the form of the field 
variables near the origin to Section [3TT| and instead here briefly state the results of this 
derivation for reference. 

We assume that the field variables m(r), S(r) and 0Oj(r) have regular Taylor series 
expansions about r = 0: 

m(r) = mo + ra\r + m27" 2 + 0(r 3 ); 
S{r) = S + S ir + S 2 r 2 + O(r 3 y, 

ujj(r) = uo jfl + LOj^r + ^ >2 r 2 + 0(r 3 ); (19) 

where the rrii, Si and oujj are constants. The constant Sq is no n- zero in order for the 
metric to be regular at the origin, but otherwise arbitrary since the field equations 
involve only derivatives of S. 

Regularity of the metric and curvature at the origin immediately gives: 

m = mi = m 2 = 0, S% = 0, uj^\ = (20) 

and 

u jt0 = ± v /j(N-j). (21) 

Without loss of generality (due to (fTTj) ). we take the positive square root in (|2TI) . 

The expansions (EEHJ) are substituted into the field equations ( JT21T151) to determine 
the values of the remaining constants in the expansions. The details of this analysis 
are presented in Section I3TTI but the upshot is as follows. We define a vector u>(r) = 
(u>i(r), uj 2 (r), . . . , uiN-i(r)) T . It turns out that, in order to have A^ — 1 independent 
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parameters to describe the N — 1 gauge field degrees of freedom O0j(r), it is necessary 
to expand the Uj(r) to order 0(r N ). Then the expansions f|T9|) become 

m(r) = m^r 3 + 0(r ); 
S{r) = S + O(r 2 ); 



N 
k=2 



^ + J2 z k( r ) v kr h + O(r N+1 ), (22) 



where 



u> = (VN=T, ^2{N -2),...,y/(N- 1)) T . (23) 

The functions Zk(r), k = 2, . . . , N are determined completely by the arbitrary constants 
Zk — 2fc(0) for k = 2, . . . , N. Together with A, the constants z 2 , ■ ■ ■ , determine the 
constant (see equation ( l56l) ). The A^— 1 vectors Vk are constant and their components 
are fixed (see equation ( |62i) ). The expansions ( 1221) therefore give the field variables in 
terms of the A^ + 1 parameters z 2 , . . . , z N , A and So. 

2.2.2. Event horizon For black hole solutions, we assume that there is a regular, non- 
extremal event horizon at r = rn, where /x(r) has a single zero. This fixes the value of 
m(r h ) to be 

Ar^ 

2m{r h ) = Th --±. (24) 

We assume that the field variables Wj(r), m(r) and S(r) have regular Taylor series 
expansions about r = r h : 

m(r) = m(r h ) + m'(r h ) (r - r h ) + O (r - r h ) 2 ; 
w i( r ) = ^j^h) + Uj(r h ) (r - r h ) + O (r - r h ) 2 ; 

S(r) = S(r h ) + S'(r h ) (r - r h ) + O (r - r h ) . (25) 

Setting /x(r/j) = in the Yang-Mills equations f|T2|) fixes the derivatives of the gauge 
field functions at the horizon: 

^ - -TTTTvffw (26) 

2m{r h ) - 2rfj) e (r h ) ^ 

Therefore the expansions (|25l) are determined by the A^ + 1 quantities u>j(rh), rh, S(rh) 
for fixed cosmological constant A. For the event horizon to be non-extremal, it must be 
the case that 

2m'(r h ) = 2r\ Ve {r h ) < 1 - Ar 2 h , (27) 

which weakly constrains the possible values of the gauge field functions ujj(r h ) at the 
event horizon. Since the field equations ( fT2lfT5l) are invariant under the transformation 
ffT7|) . we may consider ujj{r h ) > without loss of generality. 
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2.2.3. Infinity At infinity, we require that the metric (jl]) approaches adS, and therefore 
the field variables u>j(r), m(r) and S(r) converge to constant values as r — ► oo. We 
assume that the field variables have regular Taylor series expansions in r _1 near infinity: 

M + 0(r- 1 )- S{r) = 1 + {r- 1 ); 

+ O (r- 1 ) . (28) 



m(r) 



If the space-time is asymptotically flat, with A = 0, then the values of u)j )0O are 
constrained to be 



±y/j(N-j). (29) 



This condition means that the asymptotically flat black holes have no magnetic charge 
at infinity, or, in other words, these solutions have no global magnetic charge. Therefore, 
at infinity, they are indistinguishable from Schwarzschild black holes. However, if the 
cosmological constant is negative, then there are no a priori constraints on the values 
of u)j jCa . In general, therefore, the adS black holes will be magnetically charged. The 
fact that the boundary conditions at infinity in the A < case are less restrictive than 
those in the A = case leads to the expectation of many more solutions when A < 
compared with A = 0. 

2.3. Embedded solutions 

The field equations (112|15p are non-linear and coupled, but they do have two analytic, 
trivial solutions. 

Schwarzschild-adS Setting 

u 3 {r) = ±y/j(N-j) (30) 

for all j gives the Schwarzschild-adS black hole with 

m(r) = M = constant (31) 

We note that, by setting M = 0, pure adS space is also a solution. 
Reissner-Nordstrom-adS Setting 

Uj(r) = (32) 

for all j gives the Reissner-Nordstrom-adS black hole with metric function 

2M Q 2 Ar 2 
r r 2 3 
where the magnetic charge Q is fixed by 

Q 2 = In (N + 1) (N - 1) . (34) 
o 

Only for this value of the magnetic charge is the Reissner-Nordstrom-adS black hole 
a solution of the field equations. 



Ai(r) = l _ — + 5L - — , (33) 
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As well as these effectively abelian embedded solutions, there is also a class of 
embedded su(2) non-abelian solutions. These will turn out to play an important role in 
our proof of the existence of genuinely su(N) solutions in Section HI The construction of 
the embedded su(2) solutions in adS follows from that given in [22] in the asymptotically 
flat case. First, we write the N— 1 gauge field functions Uj(r) in terms of a single function 
u{r) as follows: 

u j (r) = ±y/j(N-j)u(r) Vj = 1, . . . , N - 1, (35) 

then define 



\ N = s l-N(N-l)(N + l). 



(36) 



Next we rescale the field variables as follows: 

R = A^V; A = A^A; m{R) = A^m(r); 

S(R) = S(r); u{R) = u(r). (37) 

Note that we rescale the cosmological constant A (this is not necessary in [22] as there 
A = 0). The field equations satisfied by m(R), S(R) and Cu(R) are then 



dm 
~dR 

SdR 


where we now have 



_ = jjlG + R 

2G 
R ] 
d 2 u 



2 ~ 



and 



G 



R?\i 



2m 
~R 



dR 2 



+ 



2m - 2R?Vft - 



2Ai? 3 



diU r 



u 2 ] Cj\ 



AR 2 



du 
dR 



1 



2R L 



(l"- 2 ) : 



(38) 



(39) 



(40) 



The equations (1581) are precisely the su(2) EYM field equations. Furthermore, the 
boundary conditions (I22|25II28I) also reduce to those for the su(2) case. This is 
straightforward to see for the boundary conditions at the horizon ( J25l) or at infinity 
(I2"H|) . At the origin, the su(2) embedded solutions are given by z 2 ^ 0, with 
z$ — Z4. = ... = Zj\[ = 0. Therefore any su(2), asymptotically adS, EYM soliton 
or black hole solution can be embedded into su(N) EYM to give an asymptotically 
adS soliton or black hole. Given that the solution space of su(2) EYM solitons and 
black holes in adS has been extensively studied 0, EED, EE] , we immediately have much 
useful information about the space of solutions of su(N) EYM soltions and black holes. 
However, our purpose in this paper is to prove the existence of genuinely su(N) (that 
is, not embedded) solutions. 
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2.4- Properties of the numerical su(N) solutions 

The numerical solution of the field equations ( fT2lTl5l) and the properties of the space 
of solutions are discussed in detail in [15] . In this section we briefly review the salient 
features of the solution space for our purposes here, and illustrate these features with 
some figures for the su(3) case. 

As in the su(2) case [91 [10], we find solutions in open sets of the parameter 
space ( [z2, ■ ■ -Zn] ( 1221) for solitons and [u>i(rh), ■ ■ ■ u)N-i{rh)] f° r black holes). For all 
negative values of the cosmological constant A < 0, we find solutions in which all the 
gauge field functions have no zeros, in a neighbourhood of the embedded pure adS 
(soliton) or Schwarzschild-adS (black hole) solution, although in practice the size of 
this neighbourhood decreases rapidly as |A| decreases. The existence of such a region 
of solutions where all gauge field functions have no zeros is proven in Section 14.31 see 
Corollary [TUJ As |A| decreases, the region of phase space where we have soliton or black 
hole solutions shrinks and breaks up, eventually becoming discrete points in the limit 
|A| -> 0. 

The main interest in this paper is solutions for which the gauge field functions have 
no zeros, and particularly the behaviour of the space of solutions for large |A|. As |A| 
increases, the region of phase space where solutions exist expands rapidly, as can be 
seen by comparing the size of the regions for A = —2 and A = —3 in Figures Q] and [2] 
(for black holes) and Figures [3] and H] (for solitons). It can also be seen in these figures 
that the regions of solutions where at least one gauge field function has at least one zero 
shrink as |A| increases. For black holes, Figure 10 in [15] shows that there are no black 
hole solutions for which at least one gauge field function has at least one zero when 
A = —5. As |A| increases still further, all the solutions which exist are such that all the 
gauge field functions have no zeros. In addition, the region of phase space where there 
are solutions continues to expand without limit as |A| increases. We will show in Section 
14.41 that for any point in the phase space, for all sufficiently large |A|, this point leads 
to a black hole or soliton solution (as applicable) for which all the gauge field functions 
have no zeros. 

We have su(2) embedded solutions fl35|) along the lines Ui(r h ) = u>2(r h ) for black 
hole solutions (indicated by a dashed line in Figures HH2]) and 23 = for soliton solutions 
(not everywhere along the lines, only where the lines are contained within the region of 
solutions). The existence of genuinely su(N) solutions in a neighbourhood of embedded 
su(2) solutions, seen in Figures [THU is shown in Section l4~3l 

3. Local existence of solutions 

The first part of our proof of the existence of regular soliton and black hole solutions 
of the su(N) EYM equations ( [121151) is showing that regular solutions exist in 
neighbourhoods of the points r = 0, r = and r = 00, where the field equations 
are singular. 
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®i( r h) 

Figure 1. Solution space for su(3) black holes with = 1 and A = —2. The numbers 
of zeros of the gauge field functions, n\ and ri2, are shown for the various regions 
of the solution space. The grey region denotes those areas of the (wi(r/j), w 2 (r? l )) 
plane for which the constraint (f27|) is satisfied, but for which we do not find solutions. 
Outside the shaded region the constraint (|2T[) is not satisfied. The dashed line is 
^x(jh) = ^2(^1), along which lie the embedded su(2) solutions. 




Figure 2. Solution space for su(3) black holes with rt = 1 and A = —3. The numbers 
of zeros of the gauge field functions, n\ and ri2, are shown for the various regions 
of the solution space. The grey region denotes those areas of the (u>i(rfc), ^(r^.)) 
plane for which the constraint (|2T|) is satisfied, but for which we do not find solutions. 
Outside the shaded region the constraint (|27[) is not satisfied. The dashed line is 
wi(n,) = uj2{rh), along which lie the embedded su(2) solutions. 
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Figure 3. Solution space for su(3) solitons with A = —2. The numbers of zeros of 
the gauge field functions, n\ and U2, are shown for the various regions of the solution 
space. 




Figure 4. Solution space for su(3) black holes with A = —3. The numbers of zeros of 
the gauge field functions, n\ and ri2, are shown for the various regions of the solution 
space. 
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In each case, our method follows that employed in [17J for asymptotically flat su(2) 
EYM, and subsequently generalized in asymptotically flat space to any compact gauge 
group in [2l] . For su(N) EYM in asymptotically flat space, the paper [22j used a 
different approach to prove local existence. 

The method we employ here makes use of the following theorem, which is taken 
from [IT], but which we restate for convenience: 

Theorem 1 Ji7[ / Consider a system of differential equations for n + m functions 
u — (ui, U2, ■ ■ ■ u n ) T and v — (vx, t>2, . . . v m ) T of the form 




= x a 'fi (x,u,v) ; 



x— 1 = -nvi + x^gi (x, u, v) ; (41) 
ax 

with constants Ti > and integers (Tj, q > 1, and let C be an open subset of M. n such 
that the functions fi, % = 1, . . . n and gi, % — 1, . . . m are analytic in a neighbourhood of 
x = 0, u = c, v = for all c G C. Then there exists an n-parameter family of solutions 
of the system f^ip such that 

Ui {x) = Cl + ) , Vi (x) = O ) , (42) 

where Ui(x), i = 1, . . .n and Vi(x), i = l,...m are defined for c 6 C, for \x\ < x (c) 
(for some Xq(c) > 0) and are analytic in x and c. 

This theorem allows one to parameterize the family of solutions near a singular point 
of a set of ordinary differential equations. The key part of the method is putting the 
field equations into the form ( )4Ti) with x = at the singular point, which requires 
a change of variables. We need to consider the three singular points r = 0, r = r^ 
and r — ► oo separately. Once we have the field equations in the relevant form, it is 
then straightforward to prove local existence of solutions of the field equations with the 
required behaviour ( I221I25II28I) . The proof of the local existence of solutions near the 
origin is by far the most complicated of the three. The fact that the local solutions are 
analytic in the parameters Cj as well as the independent variable x will turn out to be 
extremely useful in Section HI 



3.1. Local existence of solutions near the origin 

The local existence of solutions of the field equations near the origin is technically rather 
complicated. Local existence near the origin has been proved in the asymptotically flat 
case for su(iV) gauge fields in [22], using a method which is different from that employed 
here, and for general compact gauge groups in [21] , whose method we follow. As might 
be expected, the inclusion of a negative cosmological constant A does not significantly 
change the analysis. However, the analysis is highly involved and in particular it is 
important to verify exactly how the cosmological constant appears in the derivation. 
For this reason, in this section we have included sufficient detail to make it clear where 
our analysis differs from that in [2TJ [22] . 
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We begin by following [22j and define 

,/i.v-./). ./ o v. 

Then we define new variables Uj(r) and qj(r) as follows [22J: 

i 

Uj(r) = M i (r)7 2 ; j = 1, . . . , N — 1; 



(43) 



Qj{ r ) = - Uj_ x -N-l + 2j; j 



N. 



(44) 
(45) 



Next define a new independent variable x by 

r 

an 

where A a? is defined in (|36|) . We then write the cosmological constant A and metric 
function m(r) in terms of new quantities A and m as follows: 

A 



m(r) = \f}fh{x)\ 



A 



\ 2 ' 

N 



(46) 



and from now on consider Uj (r 



Uj(\]yx) and S(r) = S(Xnx) as functions of x instead 
of r. In terms of these variables the field equations become [22] 

alS 2SG 

dx 

2Ax 3 



dm 


1 - 


dx 







2 d 


= X fl- 



HG + P 



\ 2 N x' 



dx 2 



+ 



2m 



2x_ ~ 
X 2 



C?M 7 ' 1 r 

^ + 2 



qj\Uj, 



(47) 
(48) 



where 



G 



AT-l 

En 



P 



1 r 

3=1 



(49) 



The remainder of our analysis of the solutions of the field equations near the origin is 
split into two parts. Firstly, in Section I3.1.1[ we derive the power series expansion (1221) 
for the field variables near the origin. Secondly, in Section 13.1.21 we prove the local 
existence of solutions of the field equations having this form near the origin. In this 
section, we will not be using the Einstein summation convention, and all summations 
are written explicitly. 



3.1.1. Derivation of the power series expansion In this section we seek a power series 
expansion for m and Uj in a neighbourhood of the origin x = 0, in the form 



u ; 



k x k . 



(50) 



m(x) = > m k x ; 

k=0 k=0 

A straightforward analysis of the field equations shows that mo, mi, m2 and v,j t i are all 



zero and that u 



3,0 



We seek to substitute these power series expansions into the scaled field equations 
(I47f48p . First we require the power series expansions for the quantities qj, G and P in 



su(N) EYM solitons and black holes with A < 15 

terms of the m k and Uj jk . These do not depend on A and so have the same form as in 
[22]: 

oo oo oo 

qj = J2^ xk > G = Y,G k x k ; P = Y,PkX h ; (51) 



where 



N-l fc-1 

G k = ^ ^ 7j (I + 1) (fe - i + 1) u^ +1 Uj,jfe_i + i; 
i=l £=1 

j AT fc 

7 £ <lj,Mj,k-e+2- (52) 



It is also helpful to have a power series expansion for /x: 



4 



H = ^2/j, k x k , (53) 

fc=0 

where 

yu = 1, A*i = 0> A*2 = -2m 3 - — ; (54) 
and, for > 3, 

Hk = -2m k+ i. (55) 

It is clear that there is a minor difference in this expansion compared with the 
asymptotically flat case, due to the dependence of /12 on A. It is this small difference 
of which we must carefully keep track if we are to verify that the analysis of [22] holds 
when A ^ 0. 

Examining the coefficient of x k in the first equation (f4T|) . we find, for k > 2: 

k-2 

Here, P k depends on g 3j2 , . . . ,qj tk and hence on Uj t o, . . . , Uj tk ; the quantity Hi (£ = 
0, . . . , k — 2) depends on m 3 , . . . , m k _\ and A, and Gg (£ = 2, . . . , k) depends 
on Uj t 2, ■ ■ ■ , Uj t k- Therefore equation fl56l) determines m k+1 uniquely in terms of 
777,3, • • • j m k-i, Ujfi, . . . , Uj tk and A. Apart from the dependence on A, this is exactly 
the same situation as in |22j . 

Next we turn to the more complicated equation (14*81) . After a great deal of algebra, 
we obtain, from the coefficient of x k for k > 2, the following equation: 

bj, k = 2ljUj,k -k(k-l) Uj t k + ij+iUj+ifi - 7j-i%_i,fe ; (57) 
for j = 1, . . . , N — 1, where bj k is a complicated expression which we write as 

bj,k = c j:k + c j:k , (58) 



m k+1 



\%(k + l) 
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where 

fc-2 

Cj,k = ~ Uj,eVk,t] 

1=2 

^ fc-2 JV-1 fc-2 N-l 

Cj,k = — 2 ^i,j U i,k u i,k-£ ~ AijUi^Uj^-i 

1=2 i=l 1=2 i=l 

fc-2 k-l-2N-\ 

~E E Aj U 3,P U j,k-e-2~p] (59) 

£=2 p=2 i=l 

and we have defined further quantities VJ^ and Ajj as 

V kft = - (k - 2) (A; - 1) <5 £ , fc _ 2 + ^-m-/ + 2£ (£ - 2) m fc _ /+1 ; 
6 A N 

= lj [2S it j - 8 i+1 j - 5i^i tj ] ; (60) 

and 5ij is the usual Kronecker 5. Using the (N — 1) x (N — 1) matrix A whose entries 
are Aij we can write equation ( 1571) as 

N-l 

i=i 

This is exactly the same as the equation for the corresponding quantities bj^ in the 
asymptotically flat case [22]. However, our bj^ is slightly different from that in [22], as 
Vj t k contains A. We have examined the analysis in [22] very carefully, and found that 
Vj t k arises only in Lemma 3 in [22J, and, even in the proof of Lemma 3, the precise form 
of Vj t k is not important. Therefore Theorem 1 in [22] holds also in our case, and we 
restate it here for convenience: 

Theorem 2 12^ The recurrence relations l[5b)) and ( f^Ijj determine uniquely all the 
coefficients m k and Uj^ for k > N once N — 1 arbitrary parameters have been chosen, 
one for each equation with k — 2, . . . , N. 

We comment that in this section we have taken k = 2,...N, whereas in [22], the 
parameter k has values k — 1, . . . N — 1, and the quantity Uj^ multiplies x k+1 in [22] 
rather than x k here. Some care is therefore needed in comparing results in [22] with 
those in this section. 

The form of the coefficients Uj^ for j = 1, . . . , N — 1, k = 2, . . . , N is derived 
in exactly the same way as in [22]. The (N — 1) x (iV — 1) matrix A ( 1601) has 
eigenvalues k(k — l) for k = 2,3, ...,N (proved in [22]) with right-eigenvectors 
and left-eigenvectors cr^ , and we write the components of the vectors as v kj j for 
j = 1, . . . , N — 1. The components Vkj are given in terms of Hahn polynomials [22, 27J: 
N - 1 

^ = W^] sF2 (~ k + lj ~ j + ls k] 2 ' ~ N + 1; x) ' (62) 

where the 3F2 are hypergeometric functions. We find [22] 

%,2 = P2V2J, u j>3 = (3 3 v 3tj (63) 
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and, for k = 4, . . . , N, 

u j,k = u ),k + PkV k ,j- (64) 

Here the /3&, for k = 2, . . . , N, are arbitrary constant parameters, and u* k is a special 
solution of (I6TI) fixed by the requirement that 

cr T k u M = d k f3 k (65) 

for k = 2, . . . , N, where [22] 

j -T- _ (iV + fc-l)!(iV-fc)! 

d fc - <r fc u fc - (jv _ 1}! (jv _ 2)!jfe (jfe _ 1} {2k _ 1} , {bb) 

and is the vector with components Ui } k, ■ ■ ■ , UN-i,k- The u* k f l64j) are uniquely 
determined by Theorem [2J and the N — 1 arbitrary constants completely fix the 
coefficients Uj^ ( 1641) . 

It is convenient to expand the vector u yk in terms of the eigenvectors of the matrix 
A, by expanding the particular solutions u* k in terms of the Vi and collecting terms: 

N 

u k = Y,YkVi, k = 2,...,N. (67) 

1=2 

The analysis of [22J shows that, independent of the form of Vj t k (l6"0l) . 

Y k =0 for £ > k, (68) 

so that depends only on t) 2 , . . . , v^. This result is not trivial because it depends on 
the properties of the particular solutions u* k . Furthermore, from [22], the particular 
solutions u* k depend only on the components of the vectors for i = 2, . . . , k — 1, so 
that 

Y k k = (3 k . (69) 

Therefore the first N + 1 terms in the power series for the gauge field functions 
u = (iti, . . . , un-i) can be written as 

N k N 

u = u + J2Yl Y ^ xk = U ° + Y1 PkWA (70) 

k=2 1=2 k=2 

where u = (1, . . . , 1) T ; the Y k for £ < k are determined by ^ for i — 2, . . . , k — 1; and 
we have defined vector functions f3 k (x), k = 2, . . . , N, by 

N 

(3 k (x) = Y e k x £ ~ k v k = v k [h + Y k k +1 x + ... + Y k +N x N ~ k ] . (71) 
e=k 

The vectors f3 k (x) are proportional to v k for all values of x, and so are eigenvectors of 
the matrix A (160]) for every x. The constants (3 k , k = 2, . . . , iV in flTTl) are the same as 
those in ( 1641) . and are arbitrary, but the constants Y k+1 , . . . , Y k+N are determined by the 
(3 2 , ■ ■ ■ /3j\[ in a complicated way. Returning to the original variables, we find the power 
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series expansion for the gauge field functions Uj has the form (I22p . with the components 
of f3 k (x) determining the functions z k (r) in (1221) as follows: 



i 



Pkj( x ) = z k (r)v Kj \%^ j 3 , (72) 

where f3 k (x) = (p kjl , . ..p k ,N-i) T - 

The last remaining piece in the power series expansion of the field variables is to 
find the expansion for S. First let S = exp 5, so that ( 1471) becomes 

d5 2G 



dec 



(73) 



From this equation it is clear that 5 has a power series expansion of the form 

^o + ^rET^ (74) 

N k=2 

for some arbitrary constant Sq. This then gives a power series expansion for S of the 
form ([22]). 



3.1.2. Local existence To prove the local existence of solutions of the field equations 
(I12|I15I) in a neighbourhood of the origin r = 0, with the power series expansions 
( |22l7QII7ip derived in the previous sub-section, we follow the analysis of [21J in the 
asymptotically flat case, valid for any compact gauge group. As in the previous sub- 
section, the inclusion of a negative cosmological constant A does not significantly change 
the proof, but it is nonetheless important to carefully consider the details of each step. 

Proposition 3 There exists an N + 1-parameter family of local solutions of the field 
equations (HUES]] near r = 0, satisfying the boundary conditions (22}) . and analytic in 
A, z k , S and r. 

Proof We begin with the field equations in the form (1471481) . Defining a vector 
u = (u±, . . . , mat_i) T , we write the YM equation (|48l) as a vector equation: 



2 d 2 u 

^ dx 2 



n _ 2Ax 3 2x ~ 

2m - — -Ji p 



where the vector W = (Wi, . . . , Wn-i) has components 

N-l 

= fe+i - <?j-i) Uj = 2 Uj - u j A j,iU 2 i (76) 

i=l 

and the (N — 1) x (N — 1) matrix A is defined in ( 1601) . 

Our approach is to use the eigenvectors of the (N — 1) x (TV — 1) matrix A as a 
basis for our — 1 gauge field functions, instead of using the functions Ux, . . . , Wat-i- 
From the power series derivation in the previous subsection, we are motivated to write 
the vector u(x) as a sum over eigenvectors of the matrix A (cf. (1701) ): 

N 

u = u + J2( 3 k(x)x k ; (77) 

k=2 
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where u = (1, . . . , 1) T and the (3 k (x) are vector functions of x which satisfy 

A(3 k (x) = k (k - 1) (3 k (x) for all x. (78) 

It should be emphasized that we are making no assumptions about the (3 k (x) at this 
stage. In particular, we are not assuming that f3 k (x) are regular at the origin, nor that 
they have the form ( I7TT) . Since the N — 1 vectors Vk are eigenvectors of the matrix 
A corresponding to distinct eigenvalues, they constitute a basis of the space of iV — 1 
vectors. If we write f3 k (x) = (k{x)vk for some scalar function (k(x) (which will be a 
constant factor times the (k(%) functions defined below), then we are using the N — 1 
independent functions Ck( x ), k = 2, . . . , N, instead of the N — 1 independent functions 



u , 



0*0, j 



N — 1, to describe our gauge field. 



We now multiply both sides of (!76l) by the left-eigenvectors er k of the matrix A: 



2 -T 
X fiCT k 



d 2 U 



dx 1 



+ 



2m 



2Ax 3 



A 2 



iV 



,du 



1 

— ( 

2 



0. 



(79) 



Using the fact that the left- and right-eigenvectors of the matrix A are orthogonal 
we have 

and defining new functions (k by 



r fc/3fc, 



(80) 



51) 



we obtain 



x 2 fi 



+ 



x 



kd 2 (k 



2m — 



dx 2 
2Ax 3 



+ 2 kx k - 1 ^ + k(k-l)x k - 2 ( k 
dx 



A 2 



X 



+ kx k- 

dx 



+ -<n w. 



(82) 



The complicated part of the analysis lies in casting W into a suitable form. 
However, the form of W depends only on the structure of the su(N) gauge field, without 
any reference to the form of the metric or the cosmological constant. Therefore we may 
appeal to the analysis of [5T] . which is valid for any compact semi-simple gauge group. 
The following general result was proved in |21j for any compact semi-simple gauge group 
and we have explicitly verified it for su(N): 

Lemma 4 J21] 

N Z 

-^k{k-l)(3 k {x)x k + ^T t x t (83) 

k=2 1=2 

for some Z G N and vectors Tg(x) which are regular at x = 0. 



The key part of this result is that the expansion of W as a power series about x = 
starts with terms of 0(x 2 ), i.e. there are no i = or i = 1 terms in this expansion. The 
precise forms of the vectors T{ are not particularly illuminating so we do not give them 
here. 



su(N) EYM solitons and black holes with A < 



20 



Using Lemma H] we have, for k = 2, . . . , N: 



— —k{k — 1) x% + & l T ^- 



1=2 



To progress further, we need the following result on cr^Te, which is proved for general 
gauge group in [21], and which we have verified directly for gauge group su(iV): 



Lemma 5 



- T 



if£<k + l. 
Using Lemma O equation therefore simplifies to 

\a T k W = -k{k-l)x% +1 + 



55) 



(86) 



=k+l 



The importance of this result will become apparent on returning to ( 1821 . Defining a 
new variable Xk(%) by Xk(x) = > dividing throughout by x k+l ji and rearranging, 
we obtain 



.r- 



2/cXfc - k {k - 1) ar 1 ^ + fc (fc - 1) aT V'Cfc 



X 2 /i 



2m 



2Ax 3 



2* - 
A 2 



Z-jfc-1 



(87) 



Next we turn to the forms of the functions G and P (jUJ). We write f3 2 as 



(3 2 {x) = $ + xf3i(x) 



v 2 



for some constant /3q and scalar function Pi(x), where v 2 is a right-eigenvector of the 
matrix A, with eigenvalue 2. Once again, we make no assumption about the properties 
of (3i. We can write j3\ as an analytic function of ( 2 using (!5I]I . but this precise expression 
is not important here. Using (!88~j) and the properties of the right-eigenvector V2, we are 
able to write G and P as follows: 

G = -N (N 2 - 1) x 2 (3 2 + x 3 G; P = -N (N 2 - l) x 2 (3 2 + x 3 P; (89) 
3 3 

where G is an analytic function of x, Ck and Xk, and P is an analytic function of x and 
the Cfc- The differential equation satisfied by m (|4Tj) then simplifies to 



We define yet another new variable a by 
1 



fiG + P 



a 



[m — m 3 x 3 ] , 



(90) 



(91) 
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where m 3 = 2/3q, and, after some algebra, we obtain the differential equation satisfied 
by a in the form 

da 

x— = —3a + xJ- a , (92) 
dx 

where T a is an analytic function in a, x, Xk, and A. 

Returning now to the YM equation (j57j) . we next write 

^ = 1 + x 2 pi. (93) 

Again, we make no assumptions about jl: it can be written as a function of the variables 
x and a, and the constants 7713 and A. When we apply Theorem [1] at the end of this 
section, we will obtain as a corollary that jl is a regular function of x, at least in a 
neighbourhood of x = 0. The YM equation ( 1H71) then takes the form 

x— — = -2kxk + xT Xk + ^ °~k T e+k+ix e , (94) 
x e=o 

where J- Xk is an analytic function of the variables. We make one last change of variables, 
defining [21] 

Xk = Xk + ^kTk+i- (95) 

The precise form of ct^t^+i is not important, suffice to say it is a polynomial in the Q 
and does not involve any powers of x. 

Then, summarizing, the field equations can be written as: 
da 

x— = — 3a + xj- a \ 
dx 

d(k _ T 
X dx ~ X ^ kJ 

dS 

x— = xFs; 
dx 

dA , . 

= 0; (96) 

where all the T functions are analytic in x, A, a, Xk and (k, in a neighbourhood of 
a = 0, x = 0. Applying Theorem [1] we then have, in a neighbourhood of the origin 
x = 0, solutions of the form 



a(x) = 0(x); C k (x) = ( ,k + 0(x); Xk = 0{x) 
S(x) = So + Oix), (97) 

where Co,fc and S are constants, and a, (k, Xk and S are analytic in x, A, S and Co.fc- 
Returning to the original variables, we have proven the existence of solutions of the form 
( T22l) in a neighbourhood of the origin, which are analytic in r and the parameters Zk 
and Sq. 
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3.2. Local existence of solutions near the event horizon 

As might be anticipated, the inclusion of the cosmological constant A does not change 
significantly the proof of the local existence of solutions near the event horizon r = r^. 
Our choice of variables follows those in [TTJ, [20j [21]. For su(2) EYM black holes in adS, 
the following result is stated in [9]. 

Proposition 6 There exists an N + 2-parameter family of local black hole solutions of 
the field equations (I2JT5\) nearr = r^ > 0, satisfying the boundary conditions (25\) . and 
analytic in r h , A, ujj{r h ) and r. 

Proof We define a new independent variable x by x = r — rh, and new dependent 
variables as follows [T71 [20l [21] : 

P = r, A = ^, Vi = w is = — = ^. (98) 
x x J 

The field equations (I12II15P then take the form 
dp 

x— — = x; 

dx 

x— = - X + xH x + F x ; 

dx 

dipj x£j 
dx A 

dS 2x 

x— = — Go; 

dx p 

x^ =0, (99) 



where 



dx 

1 2 
F x = --A p - V; 
P P A 

H\ — — — (1 + 2G) ; 
P 

l % = - ^)' r 
2C 

'U. = - —to (ioo) 

we have defined a quantity V which is a polynomial in the ipj'. 

1 N 

V = r 4 p e = - £ [W ~ V?-i -N-1 + 2j) 2 ] (101) 

3=1 

and pe, Wj and G are given in (I13lll4|ll6p respectively. In particular, G takes the form 
1 N ^ 

G=^E^ ( 102 ) 

3=1 
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The functions Fx, F^,, H\ and H^, are therefore polynomials in 1/p, 1/A, p, A, ipj, £,j 
and A. 

The equations ( 1991) are not yet in the required form f j4TT) . We require a further 
change of variables: 

X = X — F x ; (103) 

The differential equations for these new variables are: 

x— = — X + xGx] 
dx 

x^t=-i j + xGe i , (104) 

where 

We do not write out here the full forms of Gx and G^. , which are rather lengthy. However, 
from the forms of F\ and F^. f llOOl) . using the differential equations (1991) . it is clear that 
both Gx and G^. are polynomials in 1/p, 1/A, p, A, A, £j and A. 

The equations ( 11041) are now in the form (1411) required for the application of 
Theorem [TJ We have, in a neighbourhood of x = 0, solutions of the form 

p =r h + 0(x); X = 0(x); ipj = Wj-(r h ) + 0(x); 

i j = 0(x), S = S(r h )+0(x) 1 (106) 

with p, A, and all analytic in x, r^, Uj(rh), A and S(rh). Transforming back to our 
original variables, we have proven the existence of solutions of the field equations in a 
neighbourhood of a black hole event horizon r = r^, satisfying the boundary conditions 
( 1251) . and analytic in r, r^, ^(r^), A and S(rh). 



3. 3. Local existence of solutions near infinity 

For asymptotically flat su(N) solitons and black holes, the fact that the values of the 
gauge field functions as r — > oo are fixed by (1251) makes the proof of local existence of 
solutions in a neighbourhood of infinity rather complicated. In fact, as shown in [2T1122] . 
it is necessary to expand the gauge field functions 0Jj(r) to order r~ N+1 in a manner 
similar to the analysis in Section 13.11 near r = 0. In the asymptotically adS case, the 
fact that the values of the gauge field functions at infinity are not fixed a priori makes 
the proof of the local existence of solutions much easier. The result below is stated in 
the su(2) case in [9j. 
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Proposition 7 There exists an 2N -parameter family of local solutions of the field 
equations ( fTS|[73|) near r = oo, satisfying the boundary conditions $E3\) , and analytic 
in A, Uj tCO , M and r _1 . 

Proof In this case our new independent variable is x = r _1 , and we define new 
dependent variables as follows, following [TTJ [21] : 

\ = 2m = r (l- fi — J ; tpj = Uf, & = r 2 J r (107) 

Then the field equations ( Tl~2yT5l) take the form 
d\ 

x— = xf X ] 
ax 

X dx ~ X ^ 

d£j f 
X ^ =xf *> 

dS _ 4 

dx 

x^ =0. (108) 

dx 

The functions on the right-hand-side of fll08j) are given by 

JV-l 



(A \ — 



4 " A 



j'=i 

3 " 



3x 2 3Ax 3n 1 



[WjUj - 2x£j + 3\x 2 £j - 2x 3 £ j V] ; 



A A 

2V-1 

where V = x~ A pe is a polynomial in the ipj, defined in equation fllOip . Immediately 
we have that the /'s are analytic in a neighbourhood of x = 0, so applying Theorem [T] 
gives, in a neighbourhood of x = 0, solutions of the form 

A = 2M + 0(x); = ^j,oo + O(x); ^ = c,- + O(x); 

5 = S oc + 0(x 4 ); (110) 

for constants M, Uj >00 , c? and S^, and the solutions are analytic in x, A, and the 
constants. If we fix Soo = 1 for the space-time to be asymptotically anti-de Sitter, 
then the boundary conditions (j28p are satisfied, and the field variables are analytic in 
r _1 , A, M, 0Jj )O o and Cj. We note that we find additional free parameters in this case, 
corresponding to the coefficients of terms of 0(r _1 ) in the expansion of Uj near infinity. 
These parameters will not be important in our subsequent analysis. 
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4. Existence of soliton and hairy black hole solutions 

In this section we prove two existence theorems: firstly, the existence of su(N) solutions 
in a neighbourhood of any embedded su(2) solution; and secondly the existence of 
su(N) solutions for which all the gauge field functions u>j have no zeros provided |A| is 
sufficiently large. 



4-1. Regularity for jj, > 

Our strategy in proving the existence of genuinely su(N) soliton and black hole solutions 
of the field equations is to start with a local solution near the origin (in the soliton case) 
or black hole event horizon and then extend this solution out to infinity 

The focus of this section is to show that solutions can be extended for larger values 
of r provided that /i > 0. This lemma was proved in [T7] for su(2) EYM with A = 0, 
and stated for A < in [9]. For su(N) EYM with A = 0, the corresponding lemma is 
stated in [20j . Our method of proof is exactly the same as in [17], so we only sketch the 
details. 

Lemma 8 As long as fi > all field variables are regular functions of r . 



Proof Consider an interval X = [r±, r 2 ), where r± < r 2 and r\ > in the soliton case, 
or r\ > rh in the black hole case. Assume that all field variables are regular on the 
interval X and that fi(r) > for all r G X = [r±, r 2 ]. We need to show that all other field 
variables (i.e. ujj(r), uj'^r) and S{r)) remain regular at r 2 . 
From the fact that /i(r 2 ) > we obtain 



?"2 - 



Arf 



> 2m(r 2 ). 



From (fl5l) . since r 2 pg is positive, we have 

2m(r 2 ) > 2 / //(r) uf dr. 

Jrx j=1 

Now /i(r) must have a minimum in the closed interval X, so defining 
yU min = min {ii(r) : r G X} > 0; 

we find 



r r 2 N - 1 1 

/ Y,<4*<£ 
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(111) 
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Direct integration of the second Einstein equation ( TTBI) then gives that logS and 
therefore S is finite at r 2 . 

Using the Cauchy-Schwarz inequality, we have 

N-1 N-1 r -i2 m JV-1 

Yl K'( r i) - w i( r 2)] 2 = Y 

3=1 3=1 
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Since the left-hand-side of the above equation is a sum of positive terms, and the right- 
hand-side is bounded above, we deduce that each Uj(r 2 ) is finite. 

Finally, we need to show that each ujjfo) is finite. To do this, we write the Yang- 
Mills equation (fT2"]) in the form 



CI 

Integrating this first order differential equation for /mv'- gives 

1 



"}(r a ) 



S(r 2 )id(r 2 ) 



S(ri)/j l (ri)ijj' j (ri) 



r2 SWiUi , ' 
dr 



(116) 



;ii7) 



so that oo'j(r2) is finite and we have shown that all field variables are finite on the closed 
interval X. 



4-2. Asymptotic behaviour as r — > oo 

One of the most important reasons for the abundance of solutions of su(2) EYM with 
A < compared with the asymptotically flat case is the difference in the behaviour 
of the Yang-Mills equations as r — > oo. Therefore we next turn to the asymptotic 
behaviour of the su(N) Yang-Mills equations (fT2"]) . 
As r —>■ oo, these take the form 

„3 

3 ' 3 ^ 



Ar 4 ,, 2Ar 3 , 

-J - —uj + W 3 uj 3 = 0. (118) 



As in the su(2) case [9] , these can be made autonomous by the change of variable 



1 / A / N 

-J- T (119) 



the Yang-Mills equations becoming 
d 2 u 



(h2 ■ 0. (120) 



The autonomous equations (I120p have critical points when 

1 2 1 

+ 2 

If all the ujj are non-zero at a critical point, their values are: 



Wm = ( 1 - Jj + -u 2 j+l + -uj)_ x ) Uj = 0. (121) 



J3 - ±y/j(N-j). (122) 

Linearization about the critical points reveals that the critical point at the origin is a 
centre, and all other critical points are saddles. 

We are in exactly the same situation as in the su(2) case: what is important is not 
so much the different critical point structure of the equation (11201) compared with the 
asymptotically flat case [20], but rather the nature of the variable r (11191) . In particular, 
when A < the variable r tends to zero as r — > oo, so that if we consider the solutions 
of the equations (I120p starting at some large value of r = r% and going all the way out to 
infinity, the corresponding values of r will lie in some small interval r 6 [0,Ti]. On the 
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other hand, for asymptotically flat solutions, the variable f which makes the Yang-Mills 
equations in flat space autonomous is f = logr, which will tend to infinity as r — > oo. 
Therefore, in the asymptotically flat case, any solution of the Yang-Mills equations must 
trace out the whole length of a phase path, whereas in the asymptotically adS case the 
solutions will only travel a small, finite distance along any one phase path. 

Including a small gravitational perturbation into the Yang-Mills equations (11181) 
will not change this key property of the variable r nor the finiteness of the phase paths 
taken by the field variables as r — > oo. 

4-3. Existence of su(N) solutions in a neighbourhood of embedded su(2) solutions 

We now have available all the ingredients we require to show the existence of genuinely 
su(N) EYM solitons and black holes. In this subsection we prove our first existence 
theorem, namely the existence of su(N) solutions in a neighbourhood of an embedded 
su(2) solution. This result relies heavily on the analyticity of the local solutions, proved 
in Section [31 

Proposition 9 Suppose there is an embedded su(2) solution of the field equations, with 
the gauge field functions all having k zeros. Then, all initial parameters in a sufficiently 
small neighbourhood of the initial parameters giving the embedded su(2) solution will 
give an su(N) solution of the field equations in which all gauge field functions have k 
zeros. 

Proof Suppose we have an embedded, non-trivial, su(2) solution of the field equations. 
If this is a soliton, the initial conditions at the origin will be given by (1221) with z 2 ^ 
and 23 = Zi = . . . = zn = 0. If this is a black hole, the initial conditions at the event 
horizon will be uji(r h ) = c^tV/i) = . . . = u N _i{r h ) ^ 0. From these initial conditions, 
we can integrate the field equations ( fT2lfT5l) all the way out to infinity to give a regular 
solution. Furthermore, the solution will satisfy the boundary conditions ( 1281) at infinity. 
For the rest of this section, we will assume that the magnitude of the cosmological 
constant A is fixed, and, if we are considering black hole solutions, that the radius of 
the event horizon, r h , is also fixed. We also suppose that the gauge field functions each 
have k zeros (they will all have the same number of zeros from ( l35i) . since this is an 
embedded su(2) solution). 

From the local existence theorems (Propositions [3] and [6] proved in Sections 13. II and 
13. 21 respectively) we know that there are solutions of the field equations locally near the 
origin or event horizon as applicable, for any values of the initial parameters z 2 , ■ ■ ■ z^ 
or u>i(rh), . . . c<jjv_i(r/j), and that these solutions are analytic in these initial parameters. 
For the embedded su(2) solution, it must be the case that jj,{r) > for all r e [r , oo) 
(where ro = for solitons and ro = for black holes). Therefore, by analyticity, for 
initial parameters close to the initial parameters giving the embedded su(2) solution, 
the corresponding local solutions will have /x(r) > for all r e [?"0) r 2] for some r 2 . By 
Lemma EJ these local solutions will also be regular on the interval [ro,r 2 ]. 
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Now fix some r\ 3> max{l,r }, such that |A| rf ^> 1 and such that, for the 
embedded su(2) solution, we have m(r 1 )/r 1 <C 1. By the analyticity argument above, 
provided the initial parameters Z2,---,Zn or ^lijh), ■ ■ ■ (as applicable) are in 

a sufficiently small neighbourhood of the initial parameters giving the embedded su(2) 
solution, the local su( N) solutions whose existence is guaranteed by Propositions [3] and 
[6] will be regular on the whole interval r 6 [?"o,ri]. That is, in this interval /i > for 
all r and the gauge field functions ujj will all have k zeros provided r\ is larger than the 
largest zero of the gauge field function uj of the embedded su(2) solution. Furthermore, 
at r 1; it will still be the case that m(r 1 )/r 1 <C 1 for these su(N) solutions as well as the 
nearby embedded su(2) solution. 

Since r x ^> max{l,r } and m{ri)/ri <C 1 for these su(N) solutions, we can use the 
asymptotically adS regime discussed in Section l4~2l Provided r\ is sufficiently large, the 
solutions will not move very far along their phase path as r increases from r\ to r — > oo. 
In particular, the gauge field functions will have no further zeros, and both the gauge 
field functions and their derivatives will not vary very much from their values at r = r\. 
This means that m(r)/r will continue to be very small as the solutions move along their 
phase path, and the asymptotic adS regime will continue to be valid. 

Therefore, we have shown the existence of genuinely su(N) soliton and black hole 
solutions in a neighbourhood of any embedded su(2) solution. The gauge field functions 
of these su(N) solutions will have the same number of zeros as the gauge field functions 
of the embedded su(2) solution. It should be emphasized that these are genuinely su(N) 
solutions rather than embedded su(2) solutions because the gauge field functions Uj(r) 
will not satisfy ( 1351 ). 

This result can be applied to the embedded Schwarzschild-adS or pure adS solutions 
(with initial parameters, ujj(rh) = a/ j (N — j) and Zj = respectively) to give, for any 
negative value of the cosmological constant A, genuinely su(N) solutions for which all 
the gauge field functions have no zeros. 

Corollary 10 For any A < 0, there exist su(N) black holes and solitons for which all 
the gauge field functions have no zeros. 

In practice, the size of this neighbourhood about the embedded pure adS or 
Schwarzschild-adS solutions in which we have su(N) soliton and black hole solutions 
where the gauge field functions have no zeros is negligibly small unless 

IfT 1 

|A| > KT^ = —N(N -1)(N+ 1). (123) 
o 

For the su(2) case, this can be seen in Figures 1 and 4 in [15]. The factor of arises 
from the embedding of the su(2) solutions in su(N) (see Section [231) . and means that 
the size of the neighbourhood decreases as N increases. This is illustrated in Figures 9 
and 13 in [15], where the size of the region of parameter space where we have nodeless 
black hole solutions with = 1 and A = — 1 can be compared for su(3) and su(4) EYM, 
the latter region being significantly smaller than the former. 
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4-4- Existence of solutions for sufficiently large |A| 

In the previous section we have proven the existence of black hole and soliton solutions 
of the su(N) EYM field equations, for any negative value of the cosmological constant 
A, such that all the gauge field functions have no zeros. However, the size of the region 
of parameter space where we have these solutions shrinks as N increases, for fixed |A|. 
In this section we shall consider the behaviour of the solutions for large |A|. We shall 
prove that, given fixed initial parameters at the origin or black hole event horizon, as 
applicable, for sufficiently large |A|, the corresponding solutions of the field equations are 
such that all the gauge field functions have no zeros. These solutions are of particular 
interest because we expect [14J that at least some of them will be linearly stable. 

Numerically (see Section [2.41 and [T5]). we find that for |A| sufficiently large, all 
soliton and black hole solutions of the field equations are such that all the gauge field 
functions have no zeros. This is not what we are able to prove in this section. Here we 
can show that, with fixed initial parameters, there are nodeless solutions for sufficiently 
large |A|, and it may well be that the magnitude of |A| required for the existence of 
nodeless solutions varies depending on the values of the initial parameters. We observe 
this numerically, for example, comparing Figures [1] and [2] in the su(3) case, we see that 
for uJi(rh) = 0.25, co^tj/i) — 1, the corresponding black hole solution is not nodeless 
when A = — 2 but is nodeless when A = —3, whereas for Ui(rh) = 0.75, ^2(^/1) = 1, the 
black hole is nodeless for both these values of A. 

We need to consider black hole and soliton solutions separately. 

4-4-1- Black holes In [9], for su(2) EYM black holes, it is shown that for any fixed 
value of the event horizon radius and any fixed value of the gauge field function on 
the horizon, uj(rh), then for all A such that |A| is sufficiently large, there is a black hole 
solution of the field equations such that the gauge field function uj{t) has no zeros. 

The extension of this result to su(N) EYM black holes, as might be expected, is 
not difficult so we just briefly sketch the derivation of the following result. 

Proposition 11 For fixed r^ and any fixed values of the gauge field functions at the 
event horizon, u)j{r h ), for |A| sufficiently large, there exists a black hole solution of the 
s\x(N) EYM field equations such that all the gauge field functions tOj{r) have no zeros. 

Proof Firstly, we note that for fixed and 0^(77,,), the constraint fT2~T|) for a regular 
event horizon is satisfied for all sufficiently large |A|. As in [9], it is helpful to define a 
length scale I by I 2 = —3/ A, and then new variables m and /}, which will be finite as 
|A| — > 00, £ — > as follows: 




(124) 
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= r 2 (f - ^ + r 2 ^j u'j + [2m - 2£ 2 r 3 p e + 2r 3 } J 3 + £ 2 W 3 u r (125) 
In the limit £ — > 0, these equations simplify considerably and have the unique solution 
m (r) = rh(r h ) = K 3 h ; S{r) = 0; Uj(r) = ujj(r h ). (126) 

We would like to extend these results to small, non-zero values of £ by analyticity, 
along the lines of the argument used in Section 14.31 The proof of Proposition [61 
specifically contains A, however, we may use the straightforward change of variables 
(cf. fB])): 

A = fX. (127) 
Then the equations are unchanged except that we include x^f- = instead of 
x lx = an d have the following equation for x^: 

x^ = -X + xH x + F x ; (128) 

where 

H x = -*(i + 2G); 
£ 2 2C 2 

F x =-+3p-—V. (129) 

The field equations are then all regular as £ — > 0, and the local existence result 
Proposition [6] carries over to give regular solutions in a neighbourhood of the event 
horizon. Furthermore, these solutions will be analytic in r^, tOj{rh) and I. 

Once we have the local existence of solutions near the event horizon, exactly the 
same argument as used in Section 14.31 shows that we have solutions for i sufficiently 
small. In particular, fixing some r± ^> r^, and fixing both the radius of the event 
horizon and the values of the gauge field functions on the event horizon, and varying 
just £, for £ sufficiently small we will have local solutions near the event horizon which 
are regular for all r e [rh, n] and for which all the gauge field functions have no zeros 
in this interval. Then, provided that r% is sufficiently large, we may use the asymptotic 
regime (see Section H~2l) . and find that, as r —>■ oo, these solutions remain regular and 
the gauge field functions will have no zeros. 

The results of this section are illustrated in Figures [5] and O where we have plotted 
example black hole solutions with the same values of the event horizon radius and 
values of the gauge field functions at the event horizon, Ui{r h ) and u 2 (r h ), but with 
A = -100 and A = -1000 respectively. 

The effect of increasing |A| by an order of magnitude can be seen by comparing these 
two figures. In both cases, the gauge field functions hardly vary at all from their values 
at the event horizon, and the difference between their values at the event horizon and at 
infinity reduces as |A| increases. In both figures log S is a small function, and it decreases 
as |A| increases. One interesting point is that, although we are considering very large 
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m(r) - m(r h ) 
o»i(r) 



co 2 (r) 



log S(r) 



10 2 



10 3 r 10 4 



Figure 5. Example of an su(3) black hole solution with = 1, A = —100, uj\[rh) = 2 
and u>2{fh) = 1- As r — > oo, the gauge field functions tend to the following limits: 
u>i — ► 2.0962, u)2 — * 0.9625. The function log S is not identically zero: at the horizon 
log 5 = 2.6246 x 10~ 3 . 
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Figure 6. Example of ansu(3) black hole solution with = 1, A = —1000, Wi(V/j) = 2 
and u>2(rh) = 1- As r — > oo, the gauge field functions tend to the following limits: 
u>\ — ► 2.0091, UJ2 — * 0.9964. The function logS" is not identically zero: at the horizon 
log 5= 2.3968 x 10~ 5 . 
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values of |A| in both these figures, the geometry is far from being Schwarzschild-adS. This 
is shown by the fact that the metric function m(r) is far from being constant, although 
the difference between its values on the event horizon and at infinity is decreasing 
as |A| increases, albeit slowly. Of course, Proposition [11] tells us that m(r) will be 
approximately constant only for "sufficiently large" |A|, and does not tell us how large 
"sufficiently large" is. It is clear that |A| will have to be very large indeed for m(r) to 
be approximately constant, although the gauge field functions become approximately 
constant for much smaller values of |A|. However, our main interest in this paper is 
proving the existence of solutions, even if in an extreme situation. 

Solitons Proving the existence of soliton solutions for any values of the initial 
parameters (z~2, . . . , zn) ( |22l) . for sufficiently large |A|, is much more difficult than the 
corresponding result in the black hole case. The reason for this is that, in the black 
hole case, there are two length scales: the radius of the event horizon and the length 
scale £ set by the negative cosmological constant (£ 2 = — 3/A), and we have r > r h for 
the region of space-time outside the event horizon. Therefore, with fixed, we can 
take I — > in a reasonably straightforward way because we know that r _1 is bounded. 
However, for soliton solutions, there is only one length scale, namely £, and we need to 
consider the whole range of values of r G [0, oo). We therefore have to be very much 
more careful in how we take the limit £ — > 0. 

For solitons, the result corresponding to Proposition [TT] is: 

Proposition 12 For any fixed values of parameters Z{, i = 2, . . . , N, which 
parameterize the gauge field functions near the origin, for |A| sufficiently large, there 
exists a soliton solution of the su(N) EYM field equations such that all the gauge field 
functions ooj(r) have no zeros. 

Proof Firstly, we rescale all the dimensionful quantities as follows: 

r = £x; m(r) = £rh(x). (130) 

In the limit £ — > 0, it is more convenient to write the Yang-Mills equations not in terms 
of the gauge field functions Uj f fl2l . but instead to write the gauge field functions in 
terms of new functions Cj( r ) : 

N N 

w(r) = u, + (kV k r k = w + ^ i k {x)l k x k v k . (131) 

k=2 k=2 

These functions C,k are essentially the same as the C,k in Section I3.1.2[ modulo a factor 
°f ^/j(N — j) and the rescaling of the independent variable x (compare equations ( )45j) 
and (I130p ). As shown in Section [3.1.21 the differential equations satisfied by the (k(x) 
([82"]) are rather complicated for general N. For = 3, they are given explicitly in [15J. 
However, in the limit £ — > 0, the field equations simplify considerably. The Einstein 
equations reduce to 

m = 0, 5=1, (132) 
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so that \i = 1 + x 2 , and the differential equations for the Ck{ x ) decouple: 

= x(l + x 2 )^| + 2[A; + (A; + l)x 2 ]^ 
v ' ax z L J ax 

+ xk(k + 1) 4- (133) 
This equation has a solution which is regular at x = 0: 

C*(z) oc 2 F X Q + |, ~; ^ + fc; -x 2 ^) ; (134) 

where 2 -Pi is a hypergeometric function. The constant of proportionality in fjl34j) is 
simply ^ ([22]) . It is straightforward to show, using the properties of hypergeometric 
functions, that the boundary conditions at infinity (128]) are satisfied. From (11311) . when 
£ = 0, all the gauge field functions are constant and given by luj = ±y/j(N — j). 
However, it is important to find the forms of the (k functions in the limit £ — > 0, even 
though they are multiplied by powers of £ and so do not affect the form of the gauge 
field functions in this limit. This is because they are important for small, non-zero £. 

The local existence result (Proposition [3]) of Section 13.1.21 carries over, with trivial 
amendments, using the scaled radial co-ordinate x rather than x. Therefore in a 
neighbourhood of the origin, we have local solutions of the field equations which are 
analytic in the parameters Zk and £. The same argument as used in Sections I4.3l and l4.4. II 
can then easily be used to show that there are soliton solutions for the field equations for 
sufficiently small £ (that is, sufficiently large |A|) for which all the gauge field functions 
have no zeros. 

As with the black hole solutions in the previous section, we illustrate this result 
with two example soliton solutions in Figures [7] and [HI Comparing these two figures, 
it can be seen how m(r) and log S(r) tend to zero as |A| increases, and the gauge 
field functions U12 approach a/2 everywhere. The metric function m(r) is decreasing in 
size more quickly for solitons as |A| increases, compared with the black hole case (cf. 
Figures [5] and [6]) . Our results in this section agree with [28], where it is shown that 
stable monopoles (in [28] both monopoles and dyons are considered, but here we have 
studied only purely magnetic monopole configurations) in su(2) EYM with A < are 
approximated by solutions on a pure adS background. There are non-trivial, nodeless, 
monopole solutions in pure adS (see Figure 1 in [28]) for the su(2) case, and, as in 
Section [2TB1 these can be embedded into su(N). In Proposition [12], we have shown the 
existence of nodeless su(N) solitons effectively approximated by solitons in pure adS 
(because rh ~ 0), but only in the limit of very large |A|, in which case the equations 
governing the YM degrees of freedom decouple ( 11331) . 

5. Conclusions 



In this paper we have studied the existence of four-dimensional, spherically symmetric, 
soliton and black hole solutions of su(N) EYM with a negative cosmo logical constant. 
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Figure 7. Example of an su(3) soliton solution with A = —100, z-i — —2 and Z3 = — 1. 
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Figure 8. Example of an su(3) soliton solution with A = —1000, Z2 = —2 and z$ = — 1. 



Numerical solutions for N = 3 and N = 4 were presented in |15j . but, of course, numerics 
can only find solutions for a small, finite number of values of N. The purpose in this 
paper was to prove the existence of four-dimensional, spherically symmetric, soliton 
and black hole solutions for any integer value of N. The general approach was briefly 
outlined in [14J, but in this paper we have presented all the detailed results. 

We began, in Section I2TT1 by describing the field equations for su(iV) EYM with a 
negative cosmological constant, the ansatz for the metric and gauge potential and the 
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boundary conditions at the origin, black hole event horizon (if there is one) and infinity. 
We considered only purely magnetic gauge fields, which are described by N — 1 gauge 
field functions cuj. We also described how any su(2) EYM soliton or black hole solution 
can be embedded to give an su(N) EYM solution for any value of N. The goal of the 
paper is to prove the existence of genuinely su(N) solutions, that is, solutions of the 
su(iV) EYM field equations which are not embedded su(2) solutions. 

The field equations for su(iV) EYM are singular at the origin, black hole event 
horizon and at infinity. In Section [3] we proved local existence of solutions in 
neighbourhoods of these singular points, using an approach in [T7J12T]. Our key existence 
results are in Section HI Firstly, we are able to prove the existence of genuinely su(N) 
solutions in a neighbourhood of an embedded su(2) solution, and, as a corollary, the 
existence of su(N) solutions for which all the gauge field functions have no zeros, for 
any negative value of the cosmological constant. Secondly, we show that, for any fixed 
values of the initial parameters at either the origin or the black hole event horizon, 
which fix the local solution in a neighbourhood of the relevant starting point, then for 
|A| sufficiently large, this local solution can be extended to infinity to give a soliton or 
black hole solution for which all the gauge field functions have no zeros. These solutions 
when |A| is large are of particular interest, because at least some of them are linearly 
stable under spherically symmetric perturbations. An outline of this result can be found 
in [Ej, and the details will be presented elsewhere [IB] . 

Our main result is therefore the existence of four-dimensional, spherically 
symmetric, asymptotically adS, soliton and black hole solutions of the su(N) EYM 
field equations with N — 1 gauge field degrees of freedom. For black holes, this means 
that there is no limit on the number of gauge field degrees of freedom (or 'hair') with 
which the black hole may be endowed. Of course, it is the stability of at least some of 
these black holes which is central to their importance for the no-hair conjecture, an issue 
to which we shall return (TBJ. This result raises the interesting question of whether there 
are solitons or black holes with an infinite number of gauge field degrees of freedom. 
There is some evidence [29] for black hole solutions of su(oo) EYM, but more work is 
required in this area. 

There are a number of generalizations of the model considered here which would 
merit further study. We have not considered topological black holes in adS, although 
solutions are known for the su(2) gauge group [23J. For su(2) EYM in adS, static, 
axially symmetric soliton [30] and black hole [31] solutions are known: it is likely that 
generalizations of these to su(N) EYM exist. Rotating su(2) EYM black holes in adS 
have also been found [32], and again one would expect generalizations to the larger gauge 
group to exist. However, the numerical challenges in finding such solutions cannot be 
underestimated. The stability of these solutions also remains an open question. In 
this paper we have restricted our attention to purely magnetic gauge fields. For su(2) 
EYM in adS, dyonic regular and black hole solutions exist [10], even though non-trivial 
su(2) EYM solutions in asymptotically flat space must be purely magnetic [33]. Since 
su(N) solutions in asymptotically flat space may have electric as well as magnetic charge 
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[25] . it would not be surprising to find dyonic su(N) EYM solutions in adS, although 
including the electric part of the gauge field will only introduce a single additional 
gauge field degree of freedom. The stability of the dyonic solutions, even in the su(2) 
case, remains an open question. Next, EYM solitons and black holes in more than 
four space-time dimensions have recently received much attention in the literature (see 
|34] for a review and references), and analogue higher dimensional solutions of su(N) 
EYM would be expected. Finally, one would expect to be able to extend our results to 
arbitrary compact gauge group: in [21] local existence theorems are proved for compact 
gauge group in asymptotically flat space, and extending these to include a negative 
cosmological constant should be possible, leading to existence theorems along the lines 
of those proved in Section HI 

Having shown that there is no limit to the amount of hair a black hole in adS 
can be given, the natural question arises of the consequences of these black hole 
solutions for gravitational physics. Firstly, there is the impact on the status of the 
"no-hair" conjecture: these black holes require an unlimited number of parameters to 
fully describe their configurations, and therefore are contrary to the "spirit" of the "no- 
hair" conjecture, namely that black holes are fundamentally simple objects. The fact 
that these black holes with unlimited amounts of hair exist in anti-de Sitter space may be 
significant, particularly in view of the adS/CFT (conformal field theory) correspondence 
[35] . It has been conjectured [36] that there are observables in the dual (deformed) CFT 
which are sensitive to the presence of black hole hair, and an adS / CFT interpretation 
of some stable seven-dimensional black holes with so (5) gauge fields is given in [37] 
(see also [38] for a discussion of non-abelian solutions in the context of the adS/CFT 
correspondence). Black holes with non-abelian gauge fields in supergravity have begun 
to attract attention recently [39] and understanding the supersymmetric analogues of 
the solutions we have found in this paper will be an important question to which we 
plan to return in the near future. 
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